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ABSTRACT 

We present far-UV HST/COS spectra of four nearby BL Lac Objects. BL Lac spectra are dominated 
by a smooth, power-law continuum which arises in a relativistic jet. However, the spectra are not 
necessarily featureless; weak, broad- and/or narrow-line emission is sometimes seen in high-quality 
optical spectra. We present detections of Lya emission in HST/COS spectra of Mrk421 [z = 0.030) 
and PKS 2005-489 {z = 0.071) as well as an archival HST/GHRS observation of MrkSOl (z = 0.0337). 
Archival HST/STIS observations of PKS 2155— 304 (z 0.116) show no Lya emission to a very low 
upper limit. Using the assumption that the broad-line region (BLR) clouds are symmetrically placed 
around the AGN, we use these measured Lya emission features to constrain either the relativistic P 
values for the ionizing continuum produced by the jet (in the ionization-bounded case) or the mass of 
warm gas (in the density-bounded case). While realistic P values can be obtained for all four cases, 
the values for Mrk421 and PKS 2155—304 are high enough to suggest that covering factors of broad- 
line-region clouds of ~l-2% might be required to provide consistency with earlier values of Doppler 
boosting and viewing angles suggested for this class of BL Lacs. This discrepancy also exists in the 
case of M 87, where the amount of Doppler boosting in our direction is expected to be minimal, again 
suggestive of a small covering factor of broad-line-region clouds. If, as these small covering factors 
might suggest, the assumptions of a density-bounded model could be more correct, then the observed 
Lya luminosities require that BLLac/FR I nuclei possess very little warm gas (10~* to 10~^ M0) as 
suggested by Guilbert, Fabian & McCray (1983). If these clouds are in pressure balance with a hotter 
(~ 10^ K) gas, the BLR contains too little mass to power the AGN by accretion alone. 
Subject headings: BL Lacertae objects: general, BL Lacertae objects: individual (Mrk421, 
PKS 2005-489, Mrk501, PKS 2155-304), quasars: emission lines, galaxies: nu- 
clei, ultraviolet: galaxies 



1. INTRODUCTION 

BL Lacertae Objects (BLLacs hereafter) are an ex- 
treme type of AGN in which the non-thermal contin- 
uum emission is thought to be relativistically-beamed 
and Doppler-boosted, overwhelming thermal sources of 
emission seen in most other AGN classes. These sources, 
together with the flat-spectrum radio quasars, are the 
most numerous and most luminous AGN at photon en- 
ergie s greater than a few hundred keV (e.g.. lAbdo et al.l 
[2OO9I ). While there are a few BLLacs that com- 
plete l y lack optical einission lines(e.g.. iRector &: Stockd 
120011: iSbarufatti et al.ll2006D . high signal-to -noise ratio 
(SNR) optical spectroscopy often detects weak but mod- 
erately luminous emission lines, especially in objects 
with spectral energy distributions (SEDs) that peak 

in the infr ared ("low frequency peaked BLLac s" or 

LBLs; e.g., [Sti ckel et al.' 'Scarpa fc Falomol [l997l : 
IRector fc Stocke 2001- Corbctt et al.. ,2003) like BL Lac 
itself (|Vermeulen et al lll995D . 

Less is known about the emission line re- 
gions of the high - energ y peaked BL Lacs (HBLs; 
iPadovani &: Giommil 119951 ). which have SEDs peaking 
in the UV to X-rays and extended r a dio em ission 
consistent with being iFanaroff fc RilevI ()1974D type 



1 CF Rlsl (IPerlman fc Stockd [Tool IRector fc Stock^ 
I2OOII: iGiroletti et all l2008f ). M ost HBLs have no de- 
tectable optical line emission ([Rector fc Stockd I2OOII: 
ISbarufatti et al.ll2006[) . It is possible that emission lines 
comparable to FRls (weak Hq;-|-[N II] and/or weak 
[O II] ) are present in many HBLs but are hidden at opti- 
cal wavelengths by the nonthermal A GN contin uum and 
by st arlight from the host galaxy (IBrowne"fc Ma r chU 
Il993f) . On the other hand, many HBLs possess optical 
spectra which contain the absorption lines and edges 
(e.g., Ca II HfcK break, G-band) typical of giant 
elliptical galaxies that are the usual hosts of these AGN 
tHutchi ngs fc Ncfi 1992; Wurtz, Stocke fc Yec 19961; 
Sc arpa et al.ll2000l: ISbarufatti.' Treves. fc>alomol l2005) . 

While the LBL/HBL distinction may be a result of pre- 
vious selection techniques (as suggested bv lCollinge et al.l 
120051 ). the observed properties of these two classes 
do vary; all HBLs and most LBLs have radio emis- 
sion consistent with being beamed FR Is ([Rector et al.l 
I2OOOI) while some LBLs are consistent with beamed 
emission from the more powerful Fanaroff-Riley class 2 
(FR2s) sources like the flat-spectrum radio qua sars (e.g., 
iBrinkmann et al.lll996l: IRector fc Stockd I2001D . Setting 
aside the few FR2-like LBLs, most BLLacs can be ex- 
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TABLE 1 

Observation Summary 



BLLac Object RA (J2000) Dec zagn Instrument tcxp Obs. Date 

Mrk421 11:04:27 +38:12:32 0.0300 COS/G130M 1738 2009, Dec 24 

COS/G160M 2404 2009, Dec 24 

PKS 2005-489 20:09:25 -48:49:54 0.0710 COS/G130M 2462 2009, Sept 21 

COS/G160M 1854 2009, Sept 21 

Mrk501 16:53:52 +39:45:37 0.0337 GHRS/G160M 29,196 1992, Feb 26-28 

PKS 2155-304 21:58:52 -30:13:32 0.116 STIS/E140M 14,244 1999, Nov 9 

14,244 2000, Sept 26 



plained as beamed FR 1 radio galaxies. It is unclear 
whether the differences betwe en LBL and HBL are a 
consequence of viewing angle (IPerlman fc Stockd [l993 



2001 



J_annuzi. Smith fc ElstonI 11994 iRector fc Stockd 
Nieppola ct al. 200^, peak frequenci es and luminosities 
(the so-called " Blazar sequence"; see Fossati et al." 1998'; 
Ghisellin i et a l. 1998; Ghisellini. Celotti & Costa mante 
2002; Cavalicrc & D'EUa 2002; Guct ta et al.l l2Q0l 
Maraschi et al.i .2008; .Ghisellini fc Tavecchiol I2008D . or 



a combination of both geometric and intrinsic factors 
(|Georganopoulos fc Marscheiill99 8!). 

Broad-line emission is an essential property of AGNs 
and can provide insights on a number of properties. Fun- 
damentally, the line luminosities represent a direct mea- 
sure of the amount of circumnuclear gas at parsec scales 
and the isotropic, ionizing photon environment. Direct 
measurement of the luminosity of the broad-line region 
(BLR) in BL Lacs can provide a new, direct measurement 
of the radiative power of these AGN and so can pro- 
vide an independent constraint on the amount of beam- 
ing. Current evidence suggests that the luminosity of 
any broad emission lines in BL Lacs may be linked to the 

tnchrotr on peak fr equency as originally sugges ted by 
aorgano poulos fc M arschcr (1998) and Ghisellin i et al.l 
998). If BLR photons are abundant, the cooling of the 
jet's synchrotron emission is enhanced greatly by the in- 
creased inverse- Comptonization of these photons by the 
jet. This increased cooling will decrease the peak fre- 
quency of the jet emission, resulting in a peak at IR 
energies or below (i.e., LBLs). However, without BLR 
seed photons, the jet cooling will be dominated by the 
synchrotron process and so is much less efficient. This 
results in a much higher-frequency synchrotron peak, in 
the UV or X-rays (i.e., HBLs). 

Is the lack of observed line emission in BLLacs, 
particularly HBLs, an intrinsic property of these less- 
luminous AGN, or is it a product of the photon en- 
vironment induced by the jet? The creation of sta- 
ble BLR clouds may be inhibited by soft X-ray spec- 
tra (jGuilbert. McCrav fc FabianI Il983f ) which are ob- 
served in only one of the two BLLac classes. In the 
case of LBLs, the X-ray emission is re latively weak and 
the X-ray spectrum is hard {a^ < 1; lUrrv et al.l 119961 ; 
iPadovani et al.l[2"004| ). In HBLs the X-rays are quite lu- 
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ergs s ^) and s oft (power-law 



spectral index in energy = — 1 t o —2.8; IPerlman et al.l 
[T996'. 2001 IPadovani et al.ll200lD . Thus, the soft X-ray 
continuum seen in HBL-type objects may be impeding 
their ability to create BLR clouds. This can be con- 
firmed by observational data using Lya, the most sensi- 
tive probe of BLR gas. 



In this Letter we present far-ultraviolet (FUV) spec- 
troscopy of four of the nearest HBL prototypes, in which 
weak Lya emission lines have been detected or sensitive 
upper limits can be set. In Section 2, we present new 
high SNR spectra from the Cosmic Origins Spectrograph 
(COS), newly-installed on the Hubble Space Telescope 
(HST) of Mrk421 and PKS 2005-489, an archival God- 
dard High Resolution Spectrograph (GHRS) spectrum of 
MrkSOl, and an archival Space Telescope Imaging Spec- 
trograph (STIS) spectrum of PKS 2155-304. We mea- 
sure Lya line widths and luminosities and reddening- 
corrected continuum luminosities and slopes for these 
HBLs. 

These observations provide the first simultaneous mea- 
surements of continuum and BLR line emission in HBLs. 
In Section 3, we use the observed continua to derive ion- 
izing luminosities and thereby predict the Lya line lu- 
minosities under the "nebular hypothesis" (ionization- 
bounded case). Comparison between the predicted and 
observed Lya luminosity allows an estimate of the beam- 
ing angle for the ionizing flux. Using the density- 
bounded assumption allows us to compute the amount of 
warm BLR gas assuming that the gas is optically-thin. 

2. BL LAC OBSERVATIONS AND ANALYSIS 

We present new COS FUV spectra of two BL Lac tar- 
gets (Mrk421 and PKS 2005-489) as well as archival 
observations of Mrk501 (GHRS) and PKS 2155-304 
(STIS). The four targets are summarized in Table 1. 

COS far-UV observations of Mrk421 and 
PKS 2005—489 were carried out during the first 
three months of COS science observations as part of the 
COS Guaranteed Time Observations (PID 11520, PI 
Green). Four exposures were made in each of the G130M 
(1135 < A < 1480 A) and G160M (1400 < A < 1795 A) 
medium-resolution gratings {R w 18, 000) for each 
target. Four central wavelength settings at each grat- 
ing dithered known instrumental features along the 
spectrum and provided continuous spectral coverage 
over 1135 < A < 1795 A (see IGreen et al.l I2OIOI: 
lOsterman et al.ll2010[ ). 

All COS exposures were reduced using CalCOS 
v2.11f. Flat-fielding, alignment, and coaddition of the 
processed exposures were carried out using IDL routines 
developed by the COS GTO team specifically for COS 
FUV datcQ. The details of the co addition method are 
discussed in iDanforth et al.l ()2010[) . Briefly, each expo- 

^ See |http : //casa . Colorado . edu/~danf orth/costools . html | 

for our coaddition and Hat-iielding aigorithm and additional dis- 
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Fig. 1.— COS observations of Mrk421 (top) and PKS 2005-489 (bottom) over the entire COS/G130M+G160M range. The flux has been 
corrected for reddening and are shown binned by 30 pixels (~ 4 resolution elements). Power law continua of the form -F(A) = 7o*(A/912)~°^ 
were fitted to line-free continuum regions and are shown as dashed curves. The narrow absorption features are mostly interstellar absorption 
lines; a few lines blueward of hya are intergalactic. The region around the rest-frame Lya (arrow) bounded by vertical dotted lines is 
expanded in Figure 2 below. 



sure was corrected for narrow, ~ 15%-opaque shadows 
from repellor grid wires. The local exposure time in 
these regions was reduced to give them less weight in 
the exposure-weighted final coaddition. Similarly, expo- 
sure times for data at the edges of the detectors was 
de-weighted. With four different central wavelength set- 
tings per grating, any residual instrumental artifacts 
from grid-wire shadows and detector boundaries have 
negligible effect on the final spectrum. 

Strong ISM absorption features in each exposure were 
aligned via cross-correlation and interpolated onto a com- 
mon wavelength scale. The wavelength shifts were typ- 
ically on the order of a resolution element 0.07 A, 
~ 17 km s~^) or less. The coadded flux at each wave- 
length was taken to be the exposure-weighted mean 
of flux in each exposure. To quantify the quality 
of the combined data, we identify line-free continuum 
regions at various wavelengths, smooth the data by 
the seven-pixel resolution element, and define S/N{= 
mean(flux)/stddev(flux)) ss 37 and w 20 for the Mrk421 
and PKS 2005—489 observations, respectively. The fully- 
reduced, coadded spectra of Mrk421 and PKS 2005-489 
are shown in Figure [TJ 

In addition to the two COS datasets, we analyze 
archival observations of the BLLac objects Mrk501 
(GHRS) and PKS 2155-304 (STIS). Mrk501 was ob- 
served for 29 kiloseconds with the G160M grating cover- 
ing 1222 < A < 1257 A at a FWHM of - 20 km s'^ and 
SNR« 15. PKS 2155-304 was observed for 28.5 ksec 
with the STIS/E140M grating (1140 < A < 1729 A), 
FWHM~ 7 km s-\ SNR« 27. Details of the GHRS and 



STIS datasets and reduction t echni ques can be found in 
IPenton. Stocke. fc Shulll (pOOOl) and llndebetouw fc Shulll 
(200'f, respectively. 

In order to accurately measure BL Lac continua, the 
observations were corrected for Galactic extinction. Tra- 
ditional measures of extinction are based on broad-beam 
radio observations of Galactic H I emission and may not 
take into account small-scale variations in this value. 
Since three out of four of our datasets cover the re- 
gion around Galactic Lya absorption, we can directly 
measure A^hi and calculate the extinction, E{B — V), 
along our sight line via the relationship E(B — V) = 
7Vhi/5.8 X 10^1 cm-2 (|Shull fc van Steenberd[T985l) . Un- 
certainties in fitted Galactic iVfji and the resulting errors 
in extinction correction dominate the subsequent uncer- 
tainties in continuum fitting The Mrk501 observations 
do not include 1216A, so we use the literature extinction 
value for this sight line. Obser ved spectra were c orrected 
for Galactic reddening via the iFitzpatrickl ()1999( ) param- 
eterization. 

Dereddened spectra were next blueshifted into the rest 
frame of the BLLac and line-free continuum regions 
were identified and fitted with a power law of the form 
F\ = ^0 (A/912)~"^ . The region near rest-frame Lya 
was not included in the fits. Flux at the rest-frame Ly- 

^ In Figure 1, the best-fit reddened power-laws show some sys- 
tematic deviation at both the long and short waveleght ends for 
unknown reasons-a single power law index may not be applicable, 
the COS fiux calibration by not be correct, or the far-UV extinc- 
tion correction may be deficient. These discrepencies do not affect 
our measurements of hya emission although the continuum flux at 
A < 912 A may be underestimated by ~ 10%. 
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man continuum is derived via a short extrapolation from 
the HST spectra. The Mrk501 GHRS data has such a 
small wavelength coverage that the continuum fit slope 
and LyC extrapolation are uncertain. 

Next, the Lya emission feature was measured. All 
three observed emission features were well-fit with sin- 
gle Gaussian profiles with free parameters Wccntroid, 
FWHM, and integrated intensity /(Lya). Total Lya 
luminosity (assuming isotropy) is then L(Lya) — 

/(Lya) 47r(czom/i?o)^ with i/o = 72kms-i Mpc'^ The 
significance of these features can be estimated as 



WSN^ 

VfWhm 



(1) 



where the scaling constant C is equal to the number of 
resolution elements per A. This comes to ~ 9a, ~ 15a, 
and ~ 23a for Mrk421, PKS 2005-489, and Mrk501, 
respectively. 

No Lya emission is seen in the spectrum of 
PKS 2155-304. If we assume FWHM- 3 A for the emis- 
sion line width we find a 4cr upper limit of 14^ < 44 mA or 
/(Lya) < 4 X 10^^^ erg cm~^ s""'^. Given the higher con- 
tinuum luminosity of PKS 2155-1-304, the Lya luminosity 
upper limit (L(Lya) < 1.1 x erg s'^) is comparable 
the detections in the other HBLs. 

All of the measured and derived quantities discussed 
above are listed in Table 2 for each of the four BL Lacs. 

Neither Mrk421 nor PKS 2005-489 show emission 
from the C IV AA 1548, 1550 doublet. Due to the lower 
effective area of COS at the red end of the far-UV range, 
these data are of lower S /N than those covering the Lya 
emission region. We derive upper limits on C IV emission 
for these two objects of /civ ^ 0.8 x 10"^"* erg cm~^ s~^ 
and /civ < 4.3 x 10'^^ erg cm^^ g-i f^j. Mrk421 and 
PKS 2005-489, respectively. The spectra of the other 
two BL Lac objects do not cover the rest-frame A ~ 1549 
A region. 

3. INFERENCES FOR THE BROAD-LINE 
REGIONS OF BL LACS 

With both continuum and emission-line measurements, 
we can estimate the covering fraction and beaming an- 
gle for BLLacs under the nebular hypothesis. To start 
we assume that the covering factor of the BLR clouds is 
100% and the clouds are spherically-distributed around 
the nucleus. From this simple geometry other geometries 
can be assumed, so that these calculations can be used 
as baseline values. First, we estimate the total ionizing 
flux in the BLR by extrapolating the power-law contin- 
uum to higher energies. Next, we predict the resulting 
BLR Lya line emission assuming that a) the ionizing 
flux is isotropic (clearly a poor assumption for BL Lacs) 
and, b) the broad line clouds are optically thick to ion- 
izing photons with unity covering factor (i.e., the "nebu- 
lar hypothesis" or ionization-bounded case) OR, c). the 
BLR clouds arc optically-thin individually and in toto 
(the density-bounded case). We will then interpret the 
difference between the predicted and observed Lya lu- 
minosities in the ionization-bounded case to determine 
if the required Doppler boosting is consistent with esti- 
mates from other considerations. And then we use the 
Lya line luminosities to infer BLR warm gas mass in the 
density-bounded case. 



3.1. Ionization-bounded Model 

Under Case-B recombination, the number of pho- 
tons emitted in a given line (e.g., H/3) can be related 
to the total number of ionizing photons by the ra- 
tio of their recombination coefficients (see Eq. 5.40 of 
lOsterbrock fc Ferlandl[200l . 



eff 
CtB 



h V 



dv. 



(2) 



where the effective a^^ = 3.03 x 10 cm s ^ and 

as = 2.59 X lO^^'^ cm^"^ s^^ for Case-B recombination 
at T ~ lO'' K. Lya luminosity L{Lya) is then related to 
L{Hl3) by the ratio of their respective specific intensities. 



jLya 



hvLya 

hvHfi 



as 

~efj 
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(3) 



under Case-B recombination. Substituting the power law 
fitted to the observed BL Lac spectra, Eq. 2 becomes 



/(//^)p,.o. = hvHp 



OiB 



1 

h V 



-(2-Qa) 



dv 



and the predicted Lya fiux is 



I{Lya)prc 



as 



I{Hfi\ 



(4) 



(5) 



The Lya emission predicted from the power-law con- 
tinuum fit over-predicts the observed Lya emission by 
2-4 orders of magnitude in all four cases (Table 2). The 
simplest interpretation of the Lya over-prediction is that 
the ionizing continuum seen by the BLR clouds is 2-4 
orders of magnitude less intense than what we see (i.e., 
beaming). In this case we can use the over-prediction 
factor {OFF) to calculate a required beaming factor 
{5 = [r(l — (3 cosO)\~^ , where the relativistic F and /3 
are as usually defined and is the "viewing angle" be- 
tween the outfiow axis and the observer's line-of-sight), 
the minimum required F of the outfiow and the maximum 
allowed viewing angle all assuming a unity covering fac- 
tor of BLR clouds as seen from the source. Since the bulk 
of the ionization of the BLR clouds is due to photons just 
blueward of the Lyman Limit, the OFF — 5^'^°'" , where 
the power-law sp ectral index in freq uency — {a\ ~\- 2) 
(see Table 2). See lUrrv fc Padovanil (|i995) Appendices A 
& B for equations and derivations. The exponent on 5 is 
appropriate for a spherical cloud morphology of the emit- 
ting region; a continuous, cylindrical morphology yields 
one less power on this exponent and thus more extreme 
values of F than we quote in Table 2. Also, these expres- 
sions may not be exact in these BL Lacs because we do 
not directly observe the flux and spectrum of the ioniz- 
ing radiation as a function of all off-axis angles; we are 
assuming that the dominant off-axis ionizing radiation 
source is the relativistic jet. 

On the other hand, for the large majority of viewing 
angles, the dominant UV/X-ray photon source may be 
the nuclear region of the AGN, not the jet. BLLacs and 
their likely parent population of FR 1 galaxies have little 
or no disk emission suggesting a very low accretion rate in 
a radiatively-inefficient accretion flow (Allen et al.„2006l) 
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Fig. 2. — Detailed spectra of the Lyo region in each BLLac. All data are dereddened, binned to ~ 15 — 20 km s~^, and shown over 
the range czagn = ±3000 km s~^. Power-law continuum fits are shown as dotted lines. Black arrows mark the systemic velocity of Lya 
using the absorption-line redshifts for these objects. Thick dashed curves show Gaussian Lya fits for the three detections and line centroid 
locations are marked with gray arrows. Full fit parameters are given in Table 2. All narrow absorption features are intervening interstellar 
or intergalactic absorption features unrelated to the AGN emission. 



onto t heir central black holes. As discussed in lAllen et alj 
(|2006[) . such a flow will also have a soft X-ray spectrum, 
albeit one that is thermal in shape, with a typical temper- 
ature of 0.5-1 keV. If the Lya BLR emission we observe 
in these BL Lacs is due to thermal accretion emission 
then the value for the OPF only sets lower limits on F 
and S since the off-axis jet emission must be even less 
luminous than what we have assumed for the values in 
Table 2. Either this thermal spectrum or the off-axis jet 
spectrum (which peaks in the far-UV - so ft X-ray with a 
soft X -ray spectrum, as described in e.g.. iPerlman et aD 
()2005( )') would be hostile to the formation of BLR clouds 
(see below). 

While the values of 5, minimum P, and maximum view- 
ing ang le in Table 2 are reasonable in the context o f some 
HBLs (|Hovatta et al.l [20091 : iGhisellini et al.llMot ). P = 
3 — 5 and 9 = 10° a re suggested from other considerations 
(|Urrv fc Padovanii ri995: Giroletti et al. 2006); P = 3 is 
suggested from comparisons in luminosity functions o f 
HBLs and FRl radio galaxies (|Urrv fc Padovanilll995D 
and P = 3 — 5 is suggested from radio core-to-extended 
flux ratios CPerlma n fc Stockel[l99l iRector et all 120001 : 
iGiroletti et a l. 20061) . Therefore the P values in Table 2 



for Mrk421 and PKS 2155-304 are somewhat too high 
to be consistent with earlier estimates. If we use the 
best previous values for P and 9 quoted above to deter- 
mine 6, we can then use the observed OPF to obtain a 
best estimate of the BLR cloud covering factor. In this 
case values of 1-2% are found. Such low covering factors 
are much more the expectation of the density-bounded 
model. However, it is equally possible that the values of 
P in the Mrk421 and PKS 2155-304 jets are significantly 
higher than estimates from entire samples of HBLs and 
so are consistent with an ionization bounded model. In- 
deed, by assuming that the ionization-bounded model is 
correct, the values of P in Table 2 are the first such esti- 
mates using this method. 

As a conceptual check of our methodology, we per- 
form the equivalent analysis on a representative un- 
beamed AGN with strong Lya emission. Seyfert galax- 
ies lack the relativistic jet of radio-loud AGN and the 
anisotropy of their emission is presumed to be low. Re- 
cent HST/COS observations of the nearby Sy 1.5 galaxy 
Mrk817 show strong Lya emission (/(Lya) = 8.1 x 
10~^^ erg cm~^ s~^) and a power-law continuum with 
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TABLE 2 

BL Lac Measured and Derived Quantities 



Quantity 


Mrk421 


PKS 2005-489 


Mrk501 


PKS 2155-304 


Units 




0.030 


0.071 


0.03366 


0.116 




log Af(HI) 


20.04 ± 0.03 


20.49 ±0.03 


20.04:=^ 


19.96 ±0.03 




E(B-V) 


0.019 ±0.002 


0.053 ± 0.004 


0.019:=^ 


0.016 ± 0.001 




S/Nrcs 


37 


20 


15 


27 


(see text) 




Continuum Fit: Fx 


= /o {A/912)-' 








L43±0.02 


1.07 ±0.04 


0.7 : 


1.13 ± 0.01 




lo 


27.2 ± 0.6 


5.01 ±0.22 


3.0 : 


13.4 ±0.2 


(lO-i* erg cm-2 s-^ A-^) 


X-ray Spectral Fit Parameters'" 


a 


2.10 - 2.58 


2.36 - 3.14 


1.41 - 2.18 


2.3 - 2.8 




b 


0.31 - 0.46 


0.27 ±0.09 


0.12 - 0.56 


-0.3 - 0.65 




hya Emission Feature 


•^centroid 


1251.9 


1301.1 


1256.9 


1357^= 


(A, observed) 




-60 


-190 


+70 




(km s-l) 


FWHM 


1.23 ±0.13 


4.25 ±0.24 


3.31 ±0.33 


3== 


(A) 


FWHM 


300 ± 30 


1050 ± 60 


820 ± 80 


740^^ 


(km s i) 


Equivalent Width 


-76 ± 7 


-467 ± 47 


-830 ± 83 


< 44 


(mA) 


Signficance Level 


9 


15 


23 


< 4 


(^) 


I(Lya) 


1.27 ±0.12 


2.38 ±0.11 


2.2 ±0.2 


^ 0.4 


(10 1* erg cm ^ s 1) 


L{Lya) 


2.37 ±0.22 


24.9 ±1.1 


5.2 ±0.3 


< 11 


(10*° erg s-l) 


Predicted Lyo Emission 


-'^(iS/")prodicted 


3.0 X 10-1° 


3.6 X 10"" 


~ 10-" 


1.0 X 10-1° 


(erg cm ^ s ^) 


Overprediction factor (OPF) 


2.4 X lO"' 


1.5 X 10^ 


~ 600 


> 2.6 X lO'' 




Doppler Factor {5)'^ 


16.6 


6.5 


~ 4 


> 13.6 




Minimum Required (F)'* 


8.3 


3.2 


~ 2 


> 6.8 




Maximum Viewing Angle*^ 


3.5 


8.8 


~ 16 


< 4 


(degrees) 



Literature value; not directly measured. 
^ X-ray spectral fit parameters as d efined in Perlman ct al. (200 ^): dN/dE oc _e1— l°g(^) ] Parameter ranges are taken from 
the literature of the last ten years (lAcciar i et al... 2009; ,Brinkmann et al.ll200ll. 12003. 120051: IDonato et al.ll2005l: lEdelson et al. 
[2001; Foscfiini et al. 2006; Massaro ct al. 2 00SI : IPadovani et al.ll20011:IRavasio et al.ll2004l : ISembav et al.ll2002l : ITramacere et al. 

E007a,b; Zhang ct al. 2006a, b; Zhang 200^). 
Assumed quantity; see text. 
Assuming unity covering factor of BLR clouds 



parameters /n = 1.5 x 10 erg cm ^ s ^ and a\ — 1.3 
(jWinter et all I2010D . We extrapolate the AGN con- 
tinuum via the Haardt & Madau (2005) spectrum for 
Akn 120 (another nearby Seyfert galaxy) and calculate 
the total ionizing flux as above. The predicted Lya line 
emission is 2.7 x 10-"'^^ erg cm-^ s"^, a factor of only ~ 3 
higher than observed. Given that a Lyman limit break 
is not observed in this class of objects, a covering factor 
of /blr ^ 30% is reasonable. 

A second more relevant comparison for our infer- 
ences is t he case of M87, the FRl radio galaxy 
Virgo A. iSankrit. Sembach fc CanizaresI (|l99!l ) pre- 
sented an HST/Faint Object Spectrograph (FOS) spec- 
trum of the nucleus of M87 arguing that the detected 
Lya emission had to be primarily nuclear, photion- 
ized emission by comparison to a second spectrum 
taken through an aperture displaced 0.6 arcsec nearby. 
The Lya emission of M87 is at low enough redshift 
that it is significantly affected by Galactic Lya ab- 
sorption, whose influe nce can only be estimated very 
roughly. Sankrit, Sem bach k, CanizaresI |l999) estimate 
that M87's intrinsic Lya is at least a factor of two 
larger than the 9 x lO"^"' ergs cm"^ that they mea- 
sure. If we use double their value for the Lya flux 



and take the observed continuum flux extrapolated to 
912A using a. V ^ power-law typical of the 4 BL Lacs 
presented he re, a value of OPF = 100 — 200 is ob- 
tained. The ISankrit. Sembach fc Canizared ()1999[ ) ob- 
servations are consistent with this spectral index but 
have a wavelength range too small to determine an ac- 
curate a^. This large OPF is intriguing given the ob- 
served spectral properties of M87 a nd its inferred sub- 
stant ial off-axis viewing ang le (25°; E einz fc Begelmanl 
[19971) . The minimum values of F required by the large 
jet-to-counterjet ratio for the M87 jet (3 — 5) suggest that 
the continuum is not substantially beamed in our direc- 
tion [5 ^ 1.5-2 or less if F is larger) so that the OPF 
should be quite low, 5-10 times lower than what we find. 
While other modeling of apparent superluminal motions 
within the jet (Biretta. Sparks & Macchctto 1999) sug- 
gests slightly smaller viewing angles (~ 19°), somewhat 
larger values for F are required by these motions, again 
suggesting minimal beaming in our direction. Addition- 
ally, soi ne of the nuclear L ya emission could be shock- 
heated (iDopita et al.||T997l) which would further reduce 
the amount attributable to the ionizing continuum. As 
with our more extreme BL Lac Objects the high OPF 
for M87 suggests a low covering factor of a few percent 
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for the BLR clouds. Since in the case of M87 there is 
less chance for misinterpretation due to potentially large 
beaming factors, these low covering factors are again sug- 
gestive of the density-bounded case. 

3.2. Density-bounded Model 

A second physical interpretation is that the Lya over- 
prediction is not due to beaming but to the physical 
conditions in and distribution of the BLR clouds; i.e., 
the clouds themselves are optically-thin and may also 
have small covering factor so that most of the ionizing 
photons simply escape regardless of the beaming. The 
1-2% covering factors found by applying best values to 
the ionization-bounded case may be indicative that the 
density-bounded case is more correct for this class of ob- 
jects (HBLs and FR 1 radio galaxies). Up to now we have 
assumed a large-scale-isotropic distribution of optically- 
thick clouds; however they may be optically thin and/or 
anisotropic. The BLR clouds may be more highly con- 
centrated in the e quatorial re gion of the AGN and sparser 
at the poles (e.g.. lElvisll2000t ). Regardless of their geom- 
etry, the emission line diagnostic tells us only about the 
BLR clouds irradiated by the continuum we observe and 
that we assume is coming from the jet. 

All four of our objects are observed to have steep X-ray 
spectr a. The early work of iGuilbert. McCrav fc FabianI 
(jl983l ) points out that a soft X-ray source heats 
circum-nuclear gas to a single high temperature phase 
T ^ 10® K, avoiding the creation of a stable two-phase 
medium and so inhibits the formation of large warm 
clouds in the BLR. Most AGN X-ray spectra are uni- 
formly hard with spectral index a ^ — 0.7 compared to 
the much steeper spectra seen in these four (see Table 2) 
and most other HBLs. 

Clearly the neutral gas column within the beam is low 
enough that no Lyman continuum decrement is seen. All 
four targets in this study have FUSE spectra covering the 
rest-frame Lyman limit and we see no significant change 
in continuum level on either side of this break. While it 
is not a strong constraint on the BLR neutral column, 
we can say A'hi ^ few x 10^® cm~^ for three of our four 
targets, on average within the beam. (The FUSE data 
for Mrk 501 is of poor quality and we can only set a limit 
-^Hi ^ 1-5 X 10^^ cm~^ for this object.) This may be 
due to low optical depth, small covering fraction, or a 
combination of these effects, but in any case supports a 
BL Lac BLR model which is optically-thin. 

In the optically-thin regime, we can estimate the 
amount of warm gas in the BLR by assuming that there 
is one hydrogen atom for every Lya photon emitted. The 
observed Lya luminosities for these objects (Table 2) 
are so low that this assumption yields a BLR with only 
^ 10~'^ to 10~^ Mq of warm gas. Larger amounts of 
gas could be present if the few Lya emitting clouds in 
the BLR are themselves optically thick at the Lyman 
limit so that most of the mass is not emitting any Lya. 
This cloud geometry is however at variance with sim- 
ple models of advectio n-dominated or Bondi-Hoyle ac- 
cretion flows (ADAF; iNaravan fcTl 119951 : lAUen et al.l 
[2006,). Assuming pressure equilibrium between the Lya 
emitting clouds (T ~ 10^'^ K) and a dominant coronal 
phase (T ~ 10® K) suggested by the Guilbert, Fabian 
& McCray scenario, and assuming a BLR size of a few 



parsecs, yields a mass of only '--^ 0.1 Mp,. The amounts 
inferred to power the observed AGN (jCavaliere fc D'Elial 
[200l m < 10-2 Mq yr-^) are however much larger than 
the amounts that can be supplied by the BLR gas we in- 
fer to be present in these objects. There appears to be 
so little circumnuclear gas in BL Lacs that it is unclear 
how these sources power their relativistic jets. 

It is worth considering the fact that these lines 
are not very broad at all, having FWHM < 1000 
km s-^ (|Sankrit. Sembach fc CanizaresI (|1999l ) measure a 
weak, broad component to Lya in M87 that could have 
FWHM < 3000 km s'^, ahhough the continuum is 
poorly defined in their FOS spectrum). By comparison, 
typical broad lines in Seyfert galaxies have FWHM many 
times higher, and the same is also true of the b road lines 
fou nd in L BLs ( Corbctt ct al. 2000; iScarpa fc Falom^ 
[T997t IVermeulen et al.ill995f ). In fact the line widths we 
observe, are much closer to those seen in narrow, forbid- 
den lines such as [O III]. Therefore, in addition to the 
possibilities presented above, the Lya emitting material 
could be further out in these objects, at radii > 10 pc, 
more typical of NLR gas. In this case there maybe even 
less gas available for accretion power in the BLR. 

4. CONCLUSIONS 

We have measured weak, broad Lya emission in three 
nearby high-energy peaked BLLac objects (HBLs) and 
set a sensitive upper limit for a fourth. Our Lya lu- 
minosities are ~ 10*^ erg s"^ in three objects, an or- 
der of magnitude stronger than the weak, narrow optical 
emission lines that have been observed in a few HBLs 
(log L < 40 ergs s-^: lRector et al.ll2000t iSbarufatti et~aLl 
l20Cil . To the best of our knowledge this is the first 
detection of UV line emission from this rare class of 
AGN (although Lya emission was observed in M87 by 
iSankrit. Sembach fc CanizaresI [1999). BLLacs are an 
important probe of physics at the very highest energies 
and our simultaneous measurements of both continuum 
and line emission will help constrain the models of their 
structure and kinematics. 

From the observed data, we make a simple estimate 
of the number of ionizing photons produced by the 
jet and predict the Lya line emission if this entire 
energy were radiated isotropically into optically thick 
broad-line-region clouds. The predicted Lya emission is 
2-4 orders of magnitude larger than what is observed 
and we interpret this over-prediction in two possible 
ways: 

(1) The over-prediction is a symptom of relativis- 
tic beaming angle and broad line cloud covering frac- 
tion. For our two most extreme cases (Mrk 421 
and PKS 2155— 304), the beaming angles are somewhat 
smaller than predicted via other means and we needed 
to invoke a 1-2% covering factor to be consistent with 
earlie r beaming estimates for HBLs (jUrrv fc PadovaiiH 
119951 ). However, it is equally plausible that the rela- 
tivistic r for these two sources are much higher (~ 8x) 
than the HBL population as a whole. But in the case of 
M87, for which Ly a and a weak UV continuum were de - 
tected using FOS (jSankrit. Sembach fc Canizareslll999f) . 
the over-prediction factor of OPF>100 is intriguing be- 
cause little beaming is expected in our direction. In this 
case a larger value of F does not explain the large OFF 
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value since the large viewing angle means that the contin- 
uum we observe is largely unbeamed. In this case a small 
covering factor is the most likely cause of the weakness 
of the Lya emission in M87. 

(2) Small covering factors inferred for M87 and possibly 
for Mrk 421 and PKS 2155-304 are quite similar to having 
optically-thin and/or sparse clouds. So a quite plausible 
interpretation of the paucity of Lya photons coming from 
these BLLacs is that their BLRs are density-bounded, 
not ionization-bounded. This result was predicted for 
BLLac objects by Guilbert, Fabian & McCray in 1983, 
as a consequence of their steep X-ray continua. 

Conclusion (2) can be extended in two ways. First, if 
the BLR clouds are optically thin, not all of the ionizing 
radiation is captured and converted to line emission and 
the Lya luminosities are set by the total mass in warm 
gas m the BLR, 10-4 to 10"^ Mq. Assuming pressure 
balance in the BLR between these warm clouds and a hot 
medium as in the Guilbert et al. model, requires that 
there is only ~ 0.1 Mq of gas in the BLR in toto. This is 
far too little to power the AGN that we observe in these 
objects. Therefore, other fueling mechanisms need to be 
considered for BL Lac Objects and their parent popula- 
tion, FR 1 radio galaxies. One obvious way around this 
conclusion is that, while there are very few BLR clouds 
around the nucleus in BL Lacs and FR 1 radio galaxies, 
these clouds are each very optically-thick, so that a con- 
siderable amount of mass can be hidden by being shielded 
from the ionizing continuum. However, optically-thick 
clouds seems at varian ce with an ady e ction- dominated 
accretion- flow (ADAF; iNaravan fc Y1 (jl995[) ') scenario 
for fueling these sources. 

Weak, broad Lya line emission has been de- 
tected in the nucleus of the FR 1 radio galaxy M87 
(jSankrit. S cmbach & Canizarcs 1999) with HST/FOS at 
comparable line luminosity to the Mrk 421 detection re- 
ported here. However, the low resolution (which blends 
the emission with damped Galactic Lya absorption) and 
signal-to-noise (which poorly detects the UV continuum) 
is insufficient to determine the BLR parameters for M87 
conclusively. There is also the issue that near-nuclear 
emission is thou ght to be due to sho cks like in the spec- 
tra of LINERS (jPopita et all I1997D so that the photo- 
ionized Lya could be even weaker than what is observed. 
Assuming the best available values for F and viewing 
angle in the case of M87 requires small covering fac- 
tors in the ionization-bounded solution as in the cases 
of Mrk 421 and PKS 2155-304. Thus, as with the four 
BL Lacs, a density-bounded scenario should be consid- 
ered for these low-luminosity AGN which leaves the ques- 
tion of their power source still open. The nucleus of 
M87 wiU be observed in HST Cycle 18 with COS obtain- 
ing much bettter resolution an d signal-to-noise than the 
ISankrit. Sembach fc CanizaresI (fT999( ) spectrum. 

In addition to the above, one must now wonder about 
the differences between the LBL and HBL classes of 
BLLacs. As already pointed out, a few LBLs have 
considerably more luminous broad-line emissi on than 
seen in these HBLs and other LBLs (e.g., iStickel et al.l 



19931 IScarpa fc Falomol [19971: IRector fc Stockd [200l 
Corbett et al.l 120001 : iVermeulen et al.lll995[ ). more con- 
sistent with the nebular hypothesis. In a purely viewing 
angle unification of LBLs and HBLs only a very tiny 
number of BLR clouds would be irradiated by a hard X- 
ray continuum (the LBL "region" of the celestial sphere 
as seen from the nucleus) while more BLR clouds would 
see the soft X-ray continuum characteristic of HBLs. In 
the few LBLs with luminous broad emission lines it ap- 
pears likely that much larger amounts of gas are present 
in the BLR as quasi-stable broad-line clouds. If this is the 
case the off-axis ionizing continuum must be hard enough 
to permit BLR cloud formation, at variance with the 
simple viewing angle unification scheme outlined above. 
And also, given the much larger gas mass in the BLR, 
standard accretion can be a viable power source for these 
FR2-like LBLs. These differences argue against a purely 
viewing angle based unification of all LBLs and HBLs. 

Additional work is needed to flesh out the relationship 
between the two BL Lac classes, as well as the UV line- 
emission properties of BLLacs in general. Mrk 421 and 
PKS 2005— 489 were observed early in Cycle 17 as part 
of the COS Guaranteed Time Observation (GTO) pro- 
gram. Two higher-redshift BLLacs (lES 1028+511 and 
PMNJ 1103— 2329) are scheduled for observation dur- 
ing Cycle 18 and we will analyze their Lya emission 
properties. Two additional BL Lacs at unknown red- 
shift (S5 0716-^714 and lES 1553+113) wih also be ob- 
served and we will look closely for emission features by 
which the systemic redshift can be determined. The lat- 
ter object was observed briefly in 2009 and the redshift 
was constrained via i ntervening Lya absorbers (z > 0.4; 
iDanforth efall I2O1O0 . A reobservation of M87 with 
HST/COS and STIS can be used to determine if the 
Lya emission in FR Is is consistent with expectations of 
unified schemes. Detection of both Lya and C IV emis- 
sion can test whether the gas is ionized by shocks or by 
a UV continuum. The spatial resolution along the STIS 
slit will determine how much of the observed Lya and 
continuum flux is truly nuclear. 
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